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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  NO.  l86l 

ONE-DIMENSIONAL  FLOWS  OF  AN  IMPERFECT  DIATOMIC  GAS 
By  A.  J.  Eggers,  Jr. 


SUMMARY 

With  the  assumption  that  Berthelot's  equation  of  state  accounts 
for  molecular  size  and  intermolecular  force  effects,  and  that  changes 
in  the  vibrational  heat  capacities  are  given  hy  a  Planck  term, 
expressions  are  developed  for  analyzing  one— dimensional  flows  of  a 
diatomic  gas. 

The  special  cases  of  flow  through  normal  and  oblique  shocks 
in  free  air  at  sea  level  are  investigated.  It  is  found  that  up  to 
Mach  numbers  of  10  the  pressure  ratio  across  a  normal  shock  differs 
little  from  its  ideal  gas  value;  whereas  at  Mach  numbers  above  4 
the  temperature  rise  is  considerably  below  and  hence  the  density 
rise  is  well  above  that  predicted  assuming  ideal  gas  behavior. 

The  effects  of  gaseous  imperfections  on  oblique  shock  flows  are 
studied  from  the  standpoint  of  their  influence  on  the  lift  and 
pressure  drag  of  a  flat  plate  operating  at  Mach  numbers  of  10  and 
20.  The  influence  is  found  to  be  small. 


INTRODUCTION 

A  wide  variety  of  problems  in  compressible  flow  has  been 
solved  on  the  assumption  that  air  behaves  as  an  ideal  diatomic  gas. 
This  assumption  is  justified,  provided  the  pressure  and  temperature 
range  of  interest  is  small  and  near  atmospheric.  It  is  an  experi¬ 
mental  fact,  however,  that  when  air  is  subjected  to  large  changes 
in  state  at  pressures  and  temperatures  far  removed  from  atmospheric, 
it  ceases  to  obey  the  simple  gas  law  and  exhibits  other  properties 
not  characteristic  of  an  ideal  gas.  Consequently,  flow  processes 
in  which  air  is  subjected  to  these  extreme  conditions  can  be 
expected  to  depart  from  perfect  gas  behavior.  It  is  known  that 
such  flows  will  be  encountered  by  aircraft  flying  at  high  supersonic 


2 


NACA  TN  No.  l86l 


speeds  and  in  hypersonic  wind  tunnels;  hence,  the  nature  and  extent 
of  this  departure  have  become  important  considerations  in  aerody¬ 
namics. 

Classical  theories  and  experiments  have  shown  that  three 
properties  of  a  real  gas  first  cause  it  to  exhibit  characteristics 
unlike  those  of  an  ideal  gas.  These  properties  may  be  classified 
as  thermal  and  caloric  imperfections.  Thermal  imperfections  in 
the  form  of  intermole cular  forces  and  molecular  size  effects  are 
significantly  manifest  at  low  temperatures  and  high  pressures. 

Changes  in  the  vibrational  heat  capacities  become  an  important 
caloric  imperfection  at  relatively  high  temperatures.  Circumstances 
under  which  effects  of  molecular  dissociation  and  electronic  excit¬ 
ation  become  important  (e.g.,  temperatures  appreciably  above  5000  B). 

may  be  neglected  for  the  present.  Insofar  as  gases  in  equilibrium 
are  concerned,  it  is  usually  sufficient  to  account  for  intermolec- 
ular  force  and  molecular  size  effects  with  additional  terms  in  the 
equation  of  state.  Similarly,  changes  in  the  vibrational  heat 
capacities  of  the  molecules  may  be  accounted  for  with  a  function 
of  temperature  in  the  expressions  for  the  specific  heats. 


Tsien  (reference  1)  investigated  the  effects  of  gaseous 
inperfections  on  air  flows  using  Yan  der  Waals  state  equation. 
Approximate  solutions  to  the  one- dimensional  isentropic  and  normal 
shock  equations  were  obtained,  (He  points  out  two  very  early 
papers  of  limited  scope  on  the  subject  by  A.  Busemann  and 
¥.  J.  Walker.)  The  Joule-Thorns on  effect  was  neglected  in  Tsien* s 
analysis1,  however,  thus  introducing  some  error.  Donaldson 
(reference  2)  observed  this  error  and  found  that  the  differential 
equations  of  motion  could  be  integrated  to  yield  one-dimensional 
isentropic  flow  equations,  exact  to  the  accuracy  of  Yan  der  Waals 
equation.  A  comparison  of  results  obtained  with  these  equations 
and  the  flow  equations  for  a  perfect  gas  indicated  that  appreci¬ 
able  error  would  be  made  if  the  latter  expressions  were  used  to 
predict  aerodynamic  phenomena  involving  high  temperatures  or  high 
pressures . 

A  matter  considered  by  Kantrowitz  (reference  3)  i-s 
inability  of  a  gas  composed  of  polyatomic  molecules  to  instan¬ 
taneously  adjust  its  internal  energy  to  temperature  changes  at  high 


1This  error  is  reported  to  have  been  subsequently  corrected;  however, 
the  author  is  unaware  of  the  publication. 


temperatures.  This  time  lag  in  equilibrium  partition  of  energy 
occurs  primarily  in  the  vibrational  energies  of  the  molecules 
which,  unlike  the  translational  and  rotational  energies,  require 
many  collisions  per  molecule  to  become  fully  adjusted  to  a  new 
temperature.  As  pointed  out  by  Kantrowitz,  for  rapidly  acceler¬ 
ating  flows  the  resulting  heat-capacity  lag  causes  an  entropy 
increase  and  hence  a  total-head  decrease  in  a  gas.  The  net  effect 
may  be  made  negligible,  however,  by  decreasing  the  rate  of  change 
of  temperature  of  the  fluid  with  time.  In  the  special  case  of  flow 
through  a  normal  shock  wave.  Be  the  and  Teller  (reference  k)  found 
that  heat— capacity— lag  effects  caused  the  thickness  of  the  shock 
to  increase.  (The  wave  is  considered  to  extend  over  that  portion  of 
the  fluid  in  which  equilibrium  partition  of  energy  does  not  prevail.) 
As  the  Mach  number  upstream  of  the  wave  increases,  however,  the 
shock  tends  rapidly  toward  a  sharp  discontinuity.  It  may  also  be 
observed  that,  although  heat-capacity  lag  influences  flow  within 
the  shock,  it  does  not  alter  the  magnitude  of  the  entropy  increase 
across  the  shock,  as  this  is  independent  of  the  mechanism  by  which 
available  energy  is  dissipated  within  the  wave.  These  observa¬ 
tions  are  important  for  they  lend  support  to  the  simplifying 
assumptions  that  moderately  accelerating  flows  may  be  treated  as 
isentropic  except  in  the  presence  of  a  shock  wave,  in  which  case 
the  wave  may  be  treated  as  a  discontinuity  and  the  properties  on 
both  sides  calculated  as  if  there  were  no  heat-capacity  lag. 

In  the  present  paper,  one— dimensional  flow  equations  are 
obtained  with  the  aid  of  Berthelotfs  equation  of  state.  This  equa¬ 
tion,  rather  than  Van  der  Waals,  is  employed  in  order  that  somewhat 
better  estimates  of  intermolecular  force  effects  may  be  obtained. 

It  is  assumed  that  calorio  imperfections  may  be  accounted  for  with 
a  Planck  term  in  the  expressions  for  the  specific  heats.  The 
analysis  covers  both  isentropic  and  plane-shock  flows.  Within  the 
limitations  of  the  assumptions,  both  exact  and  approximate  flow 
equations  are  developed.  The  approximate  expressions  provide 
explicit  solutions  for  all  the  Important  flow  parameters,  and  are 
designed  to  utilize  data  available  on  the  theoretical  behavior  of 
an  ideal  diatomic  gas.  (See,  e.g.,  reference  5*) 


SYMBOLS 


a  local  speed  of  sound,  feet  per  second 
A  nozzle  cross-sectional  area,  square  feet 
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molecular  size  constant,  cubic  feet  per  slug 

intermolecular  force  constant,  °R,  feet  to  the  fifth  per  slug, 
second  squared 

specific  heat  at  constant  pressure,  foot-pounds  per  slug,  °R 
specific  heat  at  constant  volume,  foot-pounds  per  slug,  °R 
base  of  natural  logarithms,  2.718.  .  . 
logarithm  to  base  e 

Mach  number  (ratio  of  local  Telocity  to  local  velocity  of 
sound) 


absolute  pressure,  pounds  per  square  foot 

gas  constant  (^1715  for  j  foot-pounds  per  slug  °R 

ab  s  o  lut  e  t  empera  tur e ,  °R 


internal  energy,  foot-pounds  per  slug 


local  velocity,  feet  per  second 


specific  volume  J  ,  cubic  feet  per  slug 
external  work  performed,  foot-pounds  per  slug 
ratio  of  specific  heats  (cp/cv) 
mass  density,  slugs  per  cubic  foot 


molecular  vibrational  energy  constant,  R 


Subscripts 

stagnation  conditions 

conditions  upstream  of  shock  wave 

conditions  downstream  of  shock  wave 


ideal  gas  quantities 
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Superscripts 

*  quantities  at  the  nozzle  throat 

1  quantum-mechanical  functions 

ANALYSIS 

It  will  he  assumed  throughout  the  analysis  that  the  gas  exhibits 
no  heat-capacity  lag.  Shock— free  flows  are  considered  isentropic, 
and  flow  through  plane  shock  waves  is  assumed  to  he  adiabatic.  Justif¬ 
ication  for  the  first  two  assumptions  is  discussed  in  the  introduction. 
The  extent  to 'which  plane  shock  processes,  particularly  at  high  Mach 
numbers,  deviate  from  adiabatic  behavior  is  left  for  future  investiga¬ 
tion. 

Berthelot's  equation  of  state  is  employed  in  two  forms,  depend¬ 
ing  upon  the  desired  range  of  applicability.  In  the  development  of 
equations  for  investigating  flows  over  a  wide  range  of  Mach  numbers, 
temperatures,  and  pressures,  the  state  equation  is  used  in  its  exact 
form, 

p  _  cp^ 

"  1-bp  T  (1) 

where  b  is  the  molecular  size  constant  and  c  is  the  inter- 
molecular  force  constant.  It  will  be  noted  that  equation  (l)  differs 

from  Yan  der  Waals  equation  by  a  factor  of  i  in  the  intermolecular 

T 

force  term.  The  introduction  of  this  factor  yields  a  variation  of 
intermolecular  force  with  temperature  which  is  in  accordance  with 
experiment.  (See,  e.g.,  reference  6.)  Thus,  as  pointed  out  in 
reference  7,  close  agreement  between  theoretical  and  experimental 
values  of  compressibility  of  gases  is  obtained. 

In  the  derivation  of  approximate  flow  equations  an  expression 
of  the  form 

P  =  pET  (1+bp  -  gSg)  (2) 

is  employed.  This  equation  is,  of  course,  a  first-order  approxi¬ 
mation  to  equation  (l). 


Flow  of  a  Diatomic  Gas  Obeying  Berthelot!s  Equation  of  State 

Isentropic  flow  of  a  Berthelot  gas  may  be  studied  with  the  aid 
of  the  differential  isentropic  expansion  equation. 
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where 


and  the  Joule— Thomson  effect  is  given  b y 


(3) 


Combining  these  egressions  and  substituting  the  value  of  P  from 
equation  (l)  yields 


c„dT  -  |  dp  -  dp  .  0 


(5) 


Now  the  differential  expression  for  cy  is 


which,  upon  substituting  from  equation  (l),  may  he  integrated  to  give 

,  2cp 
Cv  -  CV  + 

where  cvf  is  a  function  which  describes  the  quantum-mechanical 
variations  of  cv  with  temperature.  The  second  term  on  the  right 
of  equation  (6)  represents  the  effects  of  gaseous  imperfections  on 

Cy  . 


The  function  chosen  for  cvf  is  determined  by  the  molecular 
structure  of  the  gas  under  consideration  and  the  temperature  range 
over  which  accurate  predictions  of  Cy*  are  desired.  For  aerody¬ 
namic  purposes ,  diatomic  gases  are  of  primary  interest.  The 
important  temperature  range  extends  from  liquefaction  temperatures 
to  several  thousand  degrees  Kankine.  A  relatively  simple  function 
may  be  written  for  cvf  in  this  case,  as  the  number  of  trans¬ 
lational  and  rotational  degrees  of  freedom  is  constant,  and  only  the 
variation  with  temperature  of  the  vibrational  heat  capacity  need  be 
considered.  This  function  is 
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(7) 


The  second  "term  in  "the  brackets,  essentially  a  Planck  term,  accounts 
for  the  vibrational  contribution  to  the  specific  heat  at  constant 
volume.  2  The  assumption  is  that  the  molecules  of  the  gas  behave  like 
linear  harmonic  oscillators  insofar  as  the  vibrational  degrees  of 
freedom  are  concerned.  (See  reference  8.) 


An  expression  governing  isentropic  expansion  of  an  imperfect 
diatomic  gas  may  now  be  obtained  by  substituting  equations  (6)  and 
(7)  Into  equation  (5)  and  integrating  from  stagnation  to  static 
states.  Since 


4 


(i )■'- 


dp  -  2pdT 

p3 


there  results  the  relation 


0  l„  (£)  +K  In  esdiiei  -  O  (£-  -  \)  +  K  In  °>-1  ^ 

1  VT  J  p(l— bp0)  TCV  fe(0/T)-l  1 


+  H 


9  e 


(0/T) 


t  (e/T)  " 
e  '  — i 


=  0 


(8) 


In  order  to  determine  the  Mach  number  of  a  stream,  it  is  neces¬ 
sary  to  find  the  velocity  of  flow  and  speed  of  sound  in  the  stream. 
These  quantities  may  be  found  by  employing  the  one-dimensional 
energy  equation. 


du  +  d(Pv)  +  VdV  = 


+  VdV  =  0 


(9) 


This  is  a  common  method  of  accounting  for  the  variation  with  temper¬ 
ature  of  the  vibrational  heat  capacities.  It  has  been  adopted  by 
Donaldson  and  others  for  imperfect  gas  studies. 
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Substituting  equations  (l),  (4),  (6),  and  (7)  into  equation  (9)  and 
integrating  from  stagnation  to  static  temperature  and  density  yields 
for  the  Telocity 


Y2  =  2. 


■jOTo-T) 


1 — © 


B0 

W7W 

—i 


(10) 


The  corresponding  speed  of  sound  is  determined  by  substituting  equa¬ 
tions  (1),  (5),  (6),  and  (7)  into  the  general  equation 


a 


2 


dP  _  /cip\  Ap'N  dT 
dp  '  \c5p/ T  +  \3f )  p  dp 


The  resulting  expression  is 


p2T 


,2  - 


ET 


2c  p 


_c_  +  R 
T2  p(l— bp) - 


(l— bp): 


'Ti 


r1+,„  _d  JYi  Y . 

e  ,  2cp\  - 

[i+(7ri,u^  i 

(0/T)]  2  ET2j  J 

1— e 

(ID 


Combining  equations  (10)  and  (ll)  yields  the  following  equation  for 
the  Mach  number. 


M2=2  l 


V 


J 


(12) 


The  specific  heats  and  the  ratio  of  specific  heats  are  readily 
obtainable  for  a  diatomic  gas  obeying  Berthelot’s  state  equation. 
The  specific  heat  at  constant  volume  cv  is  found  by  substituting 
equation  (7)  into  equation  ( 6 ),  thus  yielding 
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(13) 


The  specific  heat  at  constant  pressure  Cp  is  obtained  by  substi¬ 
tuting  equations  (l)  and  (13)  into  the  reciprocity  relation 


(dP/dT)p 
cp  =  “T  (ap/av)T 


and  setting 


The  resulting  expression  for  Cp  is 


The  ratio  of  specific  heats  7  follows  directly,  of  course,  from 
equations  (13)  and  (15)  and  may  he  written 


1  +  ^  ‘ 

/Q\  e(0/T) 

f 8  )  -  6  ^  2co! 

/2-bp  +  1  Cp  \ 

,  \1— bp  2  ET2/ 

1  1 

n 

w  iww 

w  L  -- 

2  ET2 

+  1  gcj, 

L  (l-bp)2  ET2*  - 

1  +  ( 7  -j— 1 ) 

[iff 

e(0/T)  0  1 

e  2cp 

W0/T)]2  10,2 
—  —  — 1 

(16) 


A  matter  of  special  interest  in  the  study  of  channel  flows  i8  the 
variation  of  Mach  number  with  cross-sectional  area  of  the  channel. 

The  relation  between  these  two  quantities  is  most  conveniently  deter¬ 
mined  from  the  ratio  of  cross-sectional  area  at  the  sonic  velocity 
station  (the  station  of  minimum  area  for  any  gas,  as  pointed  out  by 


10 


NACA  TN  Ho.  l86l 


Tslen)  to  croBS-sectional  area  at  a  station  of  arbitrary  Mach  number. 
The  continuity  equation  yields  this  ratio  in  the  form 


T-GO  (St)  (XT) 

With  the  aid  of  this  equation  and  equations  (l),  (8),  (ll),  and  (12), 
the  area  ratio  corresponding  to  a  particular  Mach  number  may  he 
determined  from  a  knowledge  of  the  ideal  gas  specific  heats,  character¬ 
istic  constants  of  the  gas  (i.e.,  h,  c,  and  0),  and,  for  example,  the 
stagnation  conditions-  In  general,  however,  this  computation  cannot 
he  carried  out  entirely  analytically.  For  example,  it  is  seen  from 
equation  (8)  that  none  of  the  variables  can  he  explicitly  determined; 
hence,  a  final  solution  for  any  single  variable  must  he  obtained 
graphically  or  by  equivalent  means. 

In  order  to  relate  the  fluid  properties  on  the  two  sides  of  a 
normal  shock  wave,  it  is  convenient  to  integrate  equation  (9)*  "the 
energy  equation,  across  the  wave.  The  resulting  expression  is 


Continuity  of  flow  and  conservation  of  momentum  must  also  he  satisfied 
across  the  wave.  These  requirements  may  he  expressed  analytically  in 
the  familiar  forms,  respectively. 


and 


PlVl2-P2V22=P2-Pl 


(19) 

(20) 


Equations  (l8)  through  (20 )  provide  the  additional  information 
necessary  to  determine  the  flow  through  a  normal  shock  wave.  Here 
again,  however,  final  solutions  for  certain  of  the  unknown  flow 
parameters  are  most  conveniently  obtained  graphically. 


Approximate  Flows  of  an  Imperfect  Diatomic  Gas 

This  part  of  the  study  is  concerned  with  first-order  deviations 
of  one— dimens ional  gas  flows  from  the  behavior  of  an  ideal  gas.  With 
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this  restriction,  explicit  solutions  can  he  obtained  for  all  the 
important  flow  parameters.  Equation  (2)  will,  of  course,  he  employed 
as  the  equation  of  state.  A  simplification  is  also  allowed  in  the 
temperature  function  for  the  vibrational  contribution  to  the 
specific  heats,  since  the  specific  heats  are  restricted  to  vary  only 
in  the  first  order  from  their  ideal  gas  values.  In  this  case  only 
large  values  of  9/ T  are  considered  and  the  Planck  term  may  be 
approximated  by 


Equation  (8)  may  now  be  reduced  to  the  form 


1 


This  eq.uation  can  be  readily  solved  for  p/pQ  in  terms  of  T/T0, 
retaining  only  terms  of  the  proper  order,  thus  yielding 


By  solving  for  T/T0  in  a  similar  manner,  there  is  obtained 
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A  relation  between  the  pressure  and  density  ratios  for  isentropic 
expansion  is  obtained  by  combining  equation  (2)  with  equation  (23)* 
The  resulting  expression  is 


(24) 


Equations  (22)  through  (24)  are  useful  for  determining  the  nature  and 
extent  of  small  departures  from  ideal  gas  behavior  of  a  diatomic  gas 
undergoing  isentropic  changes  in  state. 

Isentropic  relations  will  now  be  obtained  between  pressure , 
temperature  and  density  ratios,  and  the  stream  Mach  number.  Ixl  order 
to  do  this,  it  is  convenient  to  reduce  equation  (12)  to  a  form 
consistent  with  the  first-order  approximations.  Performing  this 
operation  yields 


where 

n  =  (7i2-37i+3) 


An  expression  for  the  temperature  ratio  in  ‘terms  of  Mach  number, 
stagnation  temperature  and  density,  and  gas  constants  is  obtained  "by 
substituting  equations  (2),  (22),  and  (24)  into  equation  (25)  and 
solving;  hence 


1  +  BxCMftpo+C^M)!^  +  Dx(  e-(0/To) 
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Hence,  substituting  equations  (26),  (28),  and.  (35)  into  equation  (17)  in  the  above  form  yields 
for  A* /A, 
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Substituting  equations  (2) ,  (19) ,  and  (20)  into  equation  (38) ,  dividing  through  by  equation  (35) 


and  solving  for  the  temperature  ratio  yields  the  expression: 
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In  order  to  obtain  the  relation  for  the  density  ratio  it  is  convenient  to  write 

equation  (39)  i*1  the  form 
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is  tne  stan  of  all  the  first-order  correction  terms.  The  pressure  ratio  across  a  normal 
is  now  obtained  by  combining  equations  (19) ,  (20) ,  (35)  1  and-  (44)  to  provide  the 
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Stagnation  conditions  downstream  of  the  wave  may  be  related  to  stagnation  conditions 
upstream  of  the  wave  with  the  three  identities 
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P02  (  ^2  ^  A’l/POl  ^  (49) 

P0l  \  Pi  / \PS75J  J 


in  which  all  terms  of  the  right-hand,  members  are  given  in  previously 
developed  expressions* 

It  is  important  to  note  that  the  normal  shock  expressions  are 
applicable  to  the  study  of  oblique  shock, flows.  For  example * 
equations  (39)  through  (46)  may  be  employed  directly  by  substituting 
Mi  sin  (3  for  Mx  wherever  it  occurs.  In  order  to  determine  Ms, 
however,  it  is  convenient  to  have  a  relation  between  the  wave  angle  3 
and  the  stream  deflection  angle  6.  This  relation  may  be  obtained  by 
combining  the  momentum,  energy,  and  continuity  equations  for  flow 
through  an  oblique  shock  (reference  9)  yield,  the  expression 

=  (  i 

P-L  Px  \cot  P  cot  5  +  1  / 

which  may  be  written 

cot  5  =  tan  3  7iMi2  (1— *n)  (p2/p'^)IIi  ~  1  (52) 

M2  can  then  be  calculated  from  the  expression 

M 2=([ ( Pl^  (1  -  €p  +  %)  (53) 

Mi  VTs/i  \P2A  sin  (P-6) 

which  is  obtained  from  "the  continuity  equation.  In  this  equation 
the  ideal  gas  values  for  the  temperature  and  density  ratios  are 
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those  corresponding  to  the  initial  Mach  number  Mx  and  the  wave 
angle  0.  Knowing  from  equation  (53) >  the  remaining  normal 
shock  expressions,  equations  (49)  through  (51),  can  also  be  applied 
to  oblique  shock  flows. 

It  is  observed  that,  in  general,  all  expressions  obtained  in 
the  preceding  approximate  analysis  differ  from  the  corresponding 
ideal  gas  equations  only  by  small  correction  factors  containing  the 
Mach  number,  temperature  and  density,  and  constants  of  the  gas.  As 
pointed  out  in  the  introduction,  this  simplifies  the  investigation 
of  one -dimensional  imperfect  gas  flows,  provided  effects  of  the 
imperfections  are  not  large. 


DISCUSSION  AND  CONCLUSIONS 

The  complex  manner  in  which  gaseous  imperfections  influence 
general  one-dimensional-flow  processes  precludes  the  detailed  consid¬ 
eration,  here,  of  more  than  two  flows  of  aerodynamic  interest. 

Before  discussing  these  flows,  it  may  be  worthwhile  to  establish  more 
definitely  the  conditions  under  which  imperfections  in  air  can  alter 
its  dynamical  behavior  from  that  predicted  on  the  assumption  that 
ideal  gas  laws  are  obeyed.  These  conditions  are  easily  deduced  fran 
the  equations  employed  in  the  analysis  and  a  knowledge  of  the  char¬ 
acteristic  gas  constants.  3  They  may  be  summarized  as  follows: 

1.  The  occurrence  of  temperatures  in  excess  of  800°  R,  in  which 
case  the  specific  heats  change  appreciably. 

2.  The  occurrence  of  densities  in  excess  of  10  times  sea  level 
atmospheric,  in  which  case  the  molecules  occupy  an  appreciable 
fraction  of  the  available  volume. 

3.  The  occurrence  of  sufficiently  high  pressures  and  low 
temperatures  to  cause  the  Intermolecular  forces  to  be  appreciable. 

One  additional  observation  should  be  made  concerning  the  inte¬ 
grated  effect  of  these  phenomena  on  a  particular  flow  process.  This 
is  that  the  several  imperfections  frequently  counteract  each  other 
(note,  e.g.,  the  state  equation  and  equations  (24),  (26),  and  (27)), 

3The  values  chosen  for  the  characteristic  gas  constants,  and  the 
method  by  which  they  were  determined,  are  included  in  Appendix  B. 
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and  thus  alter  certain  flow  parameters  only  slightly  from  their  ideal 
gas  values. 

The  special  cases  of  flow  through  normal  and  oblique  shocks  in 
free  air  (NACA  standard  atmospheric  conditions  at  sea  level)  will 
now  be  considered.  The  pressure  and  temperature  ratios  across  a 
normal  shock  have  been  calculated  for  Mach  numbers  from  1  to  10 
using  equations  (l),  (ll),  (l8),  (19)>  and  (20).  In  this  case  the 
terms  containing  b  and  c  are  extremely  small  (they  could,  in 
fact,  be  neglected),  and  it  is  sufficiently  accurate  to  use  the 
values  for  these  constants  given  in  Appendix  B. 

The  pressure  rise  so  determined  is  shown  in  figure  1  as  a  func¬ 
tion  of  Mach  number,  and  it  is  seen  that  it  differs  by  less  than 
6  percent  from  the  ideal  gas  value  up  to  the  maximum  Mach  number  of 
10.  This  behavior  might  be  anticipated  from  the  ideal  gas  equation 
for  the  pressure  ratio  (Appendix  A)  which  is  relatively  insensitive 
to  variations  in  7  at  high  Mach  numbers.  The  corresponding  temper¬ 
ature  rise  across  a.  normal  shock  is  illustrated  in  figure  2.  Here 
it  is  observed  that  the  effect  of  increase  in  specific  heats  appreci¬ 
ably  reduces  the  temperature  ratio  below  the  ideal  gas  values  at  the 
higher  Mach  numbers4 5.  It  may  be  deduced  from  figures  1  and  2,  that  , 
the  density  rise  will  be  considerably  above  the  ideal  gas  value  at 
the  higher  Mach  numbers.  These  results  must  be  qualified  at  Mach 
numbers  in  excess  of  7,  for  under  these  circumstances  the  temper¬ 
atures  downstream  of  the  wave  appreciably  exceed  5000°  R  which  is 
sufficiently  high  to  cause  molecular  dissociation  and  thus  possibly 
to  alter  the  flow. 

In  view  of  the  fact  that  the  changes  with  temperature  of  the 
specific  heats  of  air  do  not  strongly  influence  the  pressure  rise 
across  a  normal  shock,  up  to  Mach  numbers  of  10,  it  follows  that 
corresponding  pressure  effects  in  flow  through  oblique  shocks  will 
also  be  small,  provided  the  Mach  number  of  the  component  of  velocity 
normal  to  the  wave  is  below  10.  This  suggests,  within  the  limita¬ 
tions  observed,  that  the  effect  of  variable  specific  heats  on  those 
aerodynamic  characteristics  essentially  independent  of  viscosity  may 
be  small  for  some  shapes.  (Epstein* s  findings  in  reference  11, 
concerning  hypersonic  flows  about  bodies,  pertain  to  much  higher 
Mach  numbers  that  those  considered  here.)  The  extent  to  which  this 


4 All  ideal  gas  data  are  obtained  from  the  tables  of  reference  5. 

5 This  result  is  consistent  with  Wood’s  findings  concerning  surface 
temperatures  at  high  Mach  numbers.  (See  reference  10.) 
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is  correct  has  been  investigated  for  the  simplest  shape,  a  flat 

assuming  Prandtl-^eyer  flof about  the  leLSg  edge^  I^pr^oT' 

vm  M,  sTp  Wnd°(  Jr"  ??  e9Uati°nB  (45)  <*^*8  * 

pressure  drag  coeff icients^f "thf jLt  ^  cof  espondln8  lift  and 

numbers  of  in  on  ?V f  the  plate  were  calculated  for  Mach 
umbers  of  10  and  20,  and  for  angles  of  attack  from  0°  to  24°  and 

St,’-  T  1V°ly-  A  “'“I"1'’™  of  the  coefficients  with  those 
obtained  assuming  ideal  gas  flow  over  both  the  lower  and  utter 
surfaces  showed  that  the  effect  of  variable  specific"  Stiver, 

’  |E*S.^  the  force  coefficients  are  smaller  in  the -casp  of  +v,Q 

a“£«eo?  SLcf  )“SLbT  “Tt  than  l  P'>r0ent'  eTen  at  the  highest 
angxes  oi  attack .)  This  result  is  not  justification  for  assumW 

charnof33  tehavior  to  calculate  the  high  Mach  number  aerodynamic8 

hearL.T-1”  of,afhitrary  shapes,  other  effects  such  afZclfic 

5  tSffS  hh^cted;  but  it  does  indicate  that  the  assmpMo“ 

forces  are  obL,ined°t^  \S^?e3  °n  ”M°h  the  “Portent  pressure 
rces  are  obtained  through  oblique  or  normal  shock  compression. 


Ames  Aeronautical  Laboratory ^ 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.  ' 
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APPENDIX  A 
IDEAL  GAS  EQUATIONS 

The  following  are  ideal  gas  equations  used  in  the  analysis  of 
approximate  flows  of  an  imperfect  diatomic  gass 


Isentropic  Flow  Equations 


Normal  Shock  Equations 


2  =  (71-1)  Mx2  +  2 
1  2/iMi2  -  (ri-i) 

/Pfe\  =  27iMxg  -  (71-1) 
\Pi/ i  7i+l 
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( BZ\  -  (71+1)  Mx2 

VpJi  (7l-l)M12  +  2 

liziM  -  (7i-i)]  r(ri-i)  M  t  gi 

\T1/i  (7i+l)2  Mx2 


7i 


f  po2\  _  /^Po2\  _  r  (7i+l)  Ml2 

7i-i 

[ 27iMx2  -  (7r-l)l 

7i-i 

Vpox/  \P0l{  L^-l)  Mx2  +  2  J 

L  7i+l  J 

APPENDIX  B 

DETERMIMTION  OF  CHARACTERISTIC  GAS  CONSTANTS  FOR  AIR 

Numerical  values  for  the  characteristic  gas  constants  b  and  c 
may  be  determined  from  a  knowledge  of  the  critical  pressure  and  tem¬ 
perature  of  a  gas.  Values  determined  in  this  manner  are,  however, 
generally  less  accurate  in  a  particular  range  of  temperatures  and 
pressures  than  is  desired.  Consequently,  it  is  often  advantageous  to 
fix  b  and  c  such  that  theoretical  and  experimental  values  of  the 
primary  parameters  op,  7  and  compressibility  (i.e..  P/p)  are  in 
good  agreement  in  the  range  of  interest.  This  method  is  the  more 
suitable  of  the  two  for  most  aerodynamic  studies  as  the  variations 
in  pressure  and  temperature  are  limited;  hence,  it  has  been  employed 
in  this  paper  to  determine  b  and  c,  as  well  as  Q,  for  air. 

It  Is  sufficient  for  the  illustrative  applications  presented  here 
to  determine  the  former  constants  only  for  the  approximate  flow 
equations.  Separate  values  of  9  are  found  for  the  exact  and  approxi¬ 
mate  Planck  terms. 

As  pointed  out  previously  (note  the  development  of  approximate 
expressions  for  cp,  Cy,  and  7),  it  is  a  simple  matter  to  obtain  the 
characteristic  gas  constants  appearing  in  the  approximate  flow  equa¬ 
tions.  In  this  example,  c  and  9  were  chosen  first  to  give 
agreement  between  experimental  and  theoretical  values  of  cp  and  7, 
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the  latter  values  being  calculated  from  equations  (33)  and  (34), 
respectively.  Data  on  the  variation  of  these  quantities  with 
pressure  and  temperature  at  high  temperatures  were  obtained  from 
references  12  through  19.®  A  comparison  of  these  data  with  theory, 
setting  c  =  2.25  x  10®  °R,fts/slug,sec  and  9  =5800  R,  is  shown 
in  figures  3  and  4  for  pressures  of  0  and  144  X  10  pounds  per 
square  foot  absolute.  The  agreement  is  observed  to  "J®  8°od 
temperatures  of  3000°  R.  Using  equations  (15)  and  (l6)  to  caiculate 
cp  and  7 ,  respectively,  and  a  value  of  Q  =  5500  R,  it  is  seen  that 
at  zero  pressure  excellent  agreement  with  the  correlated  data  is 
obtained  up  to  5000°  R.  The  approximate  theory  was  then  checked  a 
low  temperatures  with  experimental  values  of  7  given  in  reference 
20.  At  temperatures  above  liquefaction  and  pressures  up  o  5 
atmospheres,  the  difference  between  theoretical  and  experimental 
values  of  7  was  found  not  to  exceed  3  percent. 

The  molecular  size  constant  b  was  chosen  to  yield  good  agree¬ 
ment  between  theoretical  and  experimental  values  of  compressibility 
of  air.  The  theoretical  values  were  calculated  from  equation  (2), 
and  the  experimental  data  were  obtained  frornTBierenoe  21.  In 
general  it  was  found  that  using  b  =  0.485  ft3/slug  and  the  previ¬ 
ously  determined  value  of  c  equation  (2)  predicts  the  compress! 
bility  of  air  to  within  2-percent  error  at  pressures  from  0  to  150 
atmospheres,  and  temperatures  from  liquefaction  to  850  R.  (The 
deviation  from  ideal  gas  behavior  is  small  above  this  temperature 
at  the  pressures  indicated. ) 
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Figure  /. —  Variation  of  pressure  ratio  across  a  normal  shock 

wave  with  Mach  number.  ( Initial  state-  NACA  standard 
atmosphere  at  sea  level.) 
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Figure  2 Variation  of  temperature  ratio  across  a  normal 
shock  wove  with  Mach  number.  (Initial  state -NACA  standard 
atmosphere  at  sea  level . ) 
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Figure  3—  Voriation  ofCp  with  temperature  at  absolute  pressures  of  0  and  144x10 
!b  /  sq  ft  ebs. 
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